Transient heat generation in a quantum dot under a step-like pulse bias 
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We study the transient heat generation in a quantum dot system driven by a step- like or a square- 
shaped pulse bias. We find that a periodically oscillating heat generation arises after adding the 
sudden bias. One particularly surprising result is that there exists a heat absorption from the 
zero-temperature phonon subsystem. Thus the phonon population in non-equilibrium can be less 
than that of the equilibrium electron-phonon system. In addition, we also ascertain the optimal 
conditions for the operation of a quantum dot with the minimum heat generation. 

PACS numbers: 65.80.-g, 44.90.-fc, 73.23.-b 



OO 



T3 

c 

o 
o 



> 

OO 
OO 

m 
o 

(N 



X 



I. INTRODUCTION 

With the technologic development in the micro- 
integration of electronic devices, the heat generation 
problem has sharply emerged. This increasing thermal 
generation greatly influences the performance and sta- 
bility of nanoscale electronic devices, and it becomes a 
main bottleneck for the future advancement of electronic 
industry. Furthermore, some familiar concepts of heat 
generation that are applicable in macroscopic systems 
may not be valid in nanoscale systems. So there is an 
imperative need to comprehend fully the heating effect 
and to try to mitigate it. 

Owing to difficulties in probing the heating process in 
nanostructures, not much progress had been made until 
several state-of-the-art experimental scenarios were con- 
ceived in recent years.— ~— There are both indirect^— and 
direct^ methods, which enable the evaluation of the ef- 
fective local temperatures of nanoscale junctions. It has 
been experimentally demonstrated that the local heat- 
ing may induce a substantial temperature increase in 
single molecular junctions^— because of the inefficient 
heat dissipation. From a microscopic point of view, the 
main factor of heat generation is the electron-phonon 
(e-p) interaction, which transfers the ordered energy of 
electric current to the disordered vibrational energy of 
atoms. In some recent measurements of differential con- 
ductance, phonon-assisted tunnelling peaks have been 
observed in various nano-devices (e.g. Cgo, carbon nan- 
otube, etc.)r~— these clearly show the existence of e- 
p interaction. As for the aspect of theory, on the one 
hand much efforts have been made to investigate the ef- 
fects of inelastic e-p scattering on electronic transport 
through different nanostructures ^r— on the other hand 
many theoretical works have devoted to understanding 
the phenomena of heat, various models and methods 
were suggested as welli^Sr— By using the nonequilibrium 
Green's function (NEGF) method, a tractable formula 
of heat generation has been derived^ But these works 
mainly focused on the steady-state cases or, at most, the 
sinusoidal-like ac bias. In order to analyze the transient 
heat generation and its dynamics, the step-like and the 
square-shaped pulse biases are considered as the opti- 
mal driving forces42r— It is not only because they can 



provide the unambiguous time scales and the transient 
behaviours of heat generation, but also because any in- 
formation in computers is presented in bytes that consist 
of a long series of "0" and "1" digits, they are realized by 
continuously switching on and off the sharp step pulse. 

In this paper, we study the transient behaviours of 
heat generation in a quantum dot (QD) system driven 
by a step-like or a square-shaped pulse bias. The QD is 
described by the Anderson model and the intra-dot e-p 
interaction is considered. By using the NEGF method, 
the transient heat generation is derived. We find a sud- 
den bias inducing a periodically oscillating heat genera- 
tion. Moreover, it can give rise to an energy transfer from 
the zero-temperature phonon subsystem to the electronic 
subsystem under certain conditions. However, after the 
bias goes back to zero and the system restores the equi- 
librium, the total heat generation at zero-temperature is 
always positive because of the third law of thermodynam- 
ics. In addition, we also discuss how to operate the QD 
with as little heat generation as possible. 

The rest of the paper is organized as follows. In Sec.lIIl 
we describe the model and deduce the formula of the 
transient heat generation. In Sec. IHII we give the numer- 
ical results and discussion. Finally, a brief conclusion is 
presented in Sec. IIVI 



II. THEORETICAL MODEL 

To begin with, we consider a single-level QD connect- 
ing to two electrode-leads and coupling to a local vibronic 
phonon mode. The whole system is described by the typ- 
ical Hamiltoniani2& 



H 



u q ala q + X q [a q + d q ]h d + e d h d 



a . k 



H.c.]. (1) 



Herein, n d — dJd, d and c a k (a = L, R) are electronic an- 
nihilation operators in the QD and lead a, e d and e a k(t) 
are the corresponding single-particle energies, and t a rep- 
resents the hopping matrix element between the QD and 
lead a. As for the phonon subsystem, a q (aj) is an annihi- 
lation (creation) operator for one specific vibronic mode 
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with frequency u q . The strength of the e-p interaction is 
described by the coupling constant X q . 

Due to the pulse bias, the electronic level t a k{t) in 
Eq. (fTJ) is time-dependent: e ak (t) — tak + A a (t), and 
A a (t) represents the external time-dependent bias po- 
tential. For the step-like pulse case, we assume that the 
system is initially in equilibrium under the zero bias; at 
t — a bias V is suddenly switched on and remains for 
t > 0, which drives the system out of equilibrium, induces 
the electric current and then generates heat. Following 
we study the transient heat generation for t > 0. Here- 
after we consider that the bias is applied solely on the left 
lead, so A R (t) = 0, A L (t) = for t < 0, and A L (t) = V 
otherwise. 

In what follows, it is rewarding to make a canonical 
transformation on the Hamiltonian H = UHU\ with 
the unitary operator U = exp[(X q /uj q )(d q — a q )d^d]: 



H = 



l k c ak + Y / t a [cl k dX + B.c.], (2) 



a,k 



where e d — ed~X 2 /uj q is the renormalized electronic level 

of the QD, and the operator X = exp[— (X q /co q )(a q — a q )]. 
We use the general formula of transient heat generation 
presented in Ref. [28l which remains unchange under the 
canonical transformation: 

Q{t) = -2Imtt g A* I dtxe- iu "^-^{ 

J — OO 

- N q {{n d (t)\n d (tiW + ((n d (t)|n d (ii)))<}. (3) 

Next we replace the operator X by its expecta- 
tion value (X) = exp{-(X q /uj q ) 2 {N q + 1/2)], which 
is reasonable under the conditions of both t a <C 
X q and t a ;§> A g i 18 ' 40 ' 41 After this approximation, 
the electron and the phonon subsystems are decou- 
pled. So we can rigorously employ the Wick's the- 
orem to deal with the two-particle Green's func- 
tion: {T^hd^hdih)]) n = {T c [h d (t)]) R {T c [h d {tx)]) R - 
(T c [d) (t)d(ti)]) R (T C [S (h)d(t)]) s . Then the heat gener- 
ation Q(t) changes into: 



Q(t) = 



2Reu q X 2 q 



- N q [G r {t, t 1 )G < (h,t) + G<(t, h)G a (h, t)} 
+ [G^t^&fat) - G < (t,t)G < (t 1 ,t 1 )}},(4) 

where G r,a,<,> (t, t\) are the standard single-particle 
Green's functions of the QD derived from the decoupled 
Hamiltonian. N q — l/[exp(u q /kBT) — l)] is the Bose dis- 
tribution function with temperature T. Then by consid- 
ering the zero temperature case (T = 0), the expression 
of transient heat generation reduces to: 

Q(t) = 2Reu q X 2 , 1 dtie - ^ 11- *) 

J —OO 

x [G < (t,t 1 )G > (t 1 ,t) -G < (t,t)G < {t 1 ,t 1 )}, (5) 



and the Green's functions G <:> take the forms asi 42 i 43 

a 

x exp[-ie(t - h) - i [ dt'A a (t')}, (6) 
G > {t,t 1 ) - G < (t,t 1 ) + G r (t,t 1 )-[G r (t,t 1 )}*. (7) 

In Eq. (J6]), the linewidth function T a — T a (X) 2 , T a = 
2irp a \t a \ 2 , with the density of state p a of lead a. /(e) = 
l/{exp[(e — fi)/kBT] + 1} is the Fermi distribution func- 
tion. To avoid trivial details, we have introduced the 
quantity A a (e,t) defined as ; 42 ' 43 



A a (e,t) = / dt'G r (t,t') 

J — OO 

xexp[ie(i-t') + i^ dt"A a {t")]. (8) 

The retarded Greenes function G r (t,t') = -i0(t - 
t') exp[— i£d(t — t') — T(t — t')] can be easily derived, in 
which r = (Tl + T-r)/2. It is independent of the bias po- 
tential A a (t). Combining the Eqs. ([3HH]), the transient 
heat generation Q(t) can be calculated for arbitrary time- 
dependent bias A a (t). Once Q(t) is known, the total heat 
generation Qr(t) from t = to the time t is straightfor- 
wardly obtained as: Q T (t) = f ^dt'Q^') = f*dt'Q(t'). 
Concretely, by taking into account of the step-like pulse 
bias, A a (e,t) is to be simplified: 

A R (e,t) = G r (e) = l/(e-? d + if), (9) 
2 L (e,t<0) = G r (e), (10) 
A L (e,t>0) = G r (e + V) + [G r {e)-G r (e + V)} 

xexp[i(e + y-e d + if)f]. (11) 

In the following calculation, we will set the chemical po- 
tential fi = of the equilibrium system at t < as the 
zero point of energy, and also assume the symmetric bar- 
riers f L = f R . 

III. RESULTS AND DISCUSSION 

Fig.QJa) shows the transient heat generation Q(t) ver- 
sus the time t with fixed V — 1 and different renormalized 
level Id- While t < 0, the bias is zero and the system is 
in equilibrium, so one always has Q(t) =0. At t = 0, 
the bias is suddenly added, then the current I(t) and the 
heat generation Q(t) emerge at t > 0. Q(t) displays a 
periodically oscillating behaviour with the period 2n/uj q . 
The oscillation is undamped and keeps all along due to 
only a discrete phonon model being considered here. Q(t) 
is positive in one half period and negative in other half 
period, which means that the transporting electron pe- 
riodically emits and absorbs phonon. Besides, the oscil- 
lating amplitude is sensitive to e^, it is the largest while 
£d = because of the occurrence of resonant tunnelling 
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FIG. 2: (Colour online) The transient behaviours in response 
to a sudden step-like pulse bias with phonon subsystem hav- 
ing a Lorentzian spectrum, (a) and (b) are Q(t) and Qr(t) 



respectively, 
and f = 1. 



The parameters are ed = 0, uj„ = 1.5, a = 0.1, 



FIG. 1: (Colour online) The transient behaviours in response 
to a sudden step- like pulse bias with phonon subsystem having 
a single vibronic mode, (a) Q(i) vs. t with V = 1. (b), (c) 
and (d) are respectively Q(t), Qr(t), and I(t) vs. t with 
id = 0. Other parameters are u) q — 1.5 and F = 1. The unit 
Jo = ef/h. 



with the largest current for that id- When id deviates 
from zero, the oscillating amplitude reduces and gradu- 
ally approaches zero. 

Next we fix the level = and investigate the be- 
haviours of Q(t) with different bias V. For the positive 
bias V, Q{t) is positive in the pre-half period and nega- 
tive in the post-half period. But for the negative V, Q(t) 
can be negative in the pre-half period [see the region of 
< t < n/uj q 0.33(27r/f) in Fig. □»]. This result 
is surprising and indicates that the suddenly-added bias 
causes the phonon to be transferred from the phonon sub- 
system to the electronic subsystem at once. Notice that 
the temperature of the system has been set equal to zero, 
so it means the phonon population can be less than the 
original zero-T value. With \V\ increasing, the oscillating 
amplitude of Q(i) also increases. For the positive-Q half 
period, the amplitude increase is very prominent. But 
for the negative-Q half period, the amplitude increase is 
small, and even decreases for the large \V\. So in the 
case of large \V\ (|V| > uj q ), more phonons are emitted 
than absorbed, which leads to a net heat generation by 
the current. 

In Fig. QJc), the total heat generation Qrif) is shown. 
Because of the oscillating Q(t), Qrif) also exhibits the 
oscillating property with time t. For the small bias 
(0 < \V\ < u) q ), Qr(t) oscillates along a horizontal line, 
in which case there is no net heat generation. But for the 
large bias (|V| > u) q ), an ascending trend will grow more 
and more noticeably, by considering the large time scale, 
Qr{t) is approximately proportional to t. In particular, 
as for — io q < V < 0, the curves of Qr{t) always lie below 
the zero point [see the curves with V = —0.5 and —1.0 
in Fig. HJc)], which means that at any time there is a net 
phonon absorption from the zero-T phonon subsystem. 



Intuitively, this seems to be impossible. But in fact, due 
to the e-p interaction, the equilibrium phonon popula- 
tion n p h = ( a \(t) a q(t)) at t < is not zero, but a finite 
value even if T = 0. By using the NEGF method and 
taking the same approximation as in the derivation of 
Q(t), 7i p h can be calculated. It is a function of Id behav- 
ing like the occupation number of electron on the QD, 
and we will mention it later in this paper (see Fig. [3]) . 
Then by choosing the same parameters in Fig. QJc), n p h 
is about 0.13A^/r 2 . From the above results, we find that 
the suddenly-applied bias can transfer a part of energy 
of such n p h phonons from the original zero-T phonon 
subsystem to the electronic subsystem. Due to the sys- 
tem at t > being in non-equilibrium, this transfer is 
not prohibited in principle. Furthermore, we have al- 
ready checked that the total heat absorption —Qr(t) 
from the phonon subsystem does not exceed the value 
of uj q n p h w 0.20A^/r in any case. 

The response of transient electric current I(t) to the 
sudden pulse bias is also shown in Fig. [IJd), m which 
I(t) quickly rises in a very short time. The rising time 
scale of I(t) is about 0.1(27r/r). However, the rise of Q(t) 
is within the time scale (1/4)(2tt /ui q ). So if we set the 
parameter T > 0.4w g , the current I(t) will emerge much 
quicker than Q{t). In this case, the device can finish 
the operation before the occurrence of heat generation. 
Additionally, in comparison with Q(t), 1) I(t) is without 
oscillation and quickly damps to a stable value because 
of the spectrum of electronic reservoirs (lead L and R) 
being continuous, 2) I(t) has the electron-hole symmetry 
[e.g. I(t,-V, r i d ) = -I(t,V,i d )}, but Q(t,-V,-e d ) f 
±Q(t,V,£d), it is because that the phonon subsystem is 
without the electron-hole symmetry. 

Up to now, we have only considered a single isolated 
phonon mode oj q in the QD. But for a real device, the 
QD usually couples to the lead and the substrate, so 
the phonon mode ui q will unavoidably acquire a small 
line width. This is to be considered in the following dis- 
cussion with the phonon subsystem having a Lorentzian 
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FIG. 3: (Colour online) The occupation number n p h of 
phonon and the occupation number n c of electron on the QD 
at zero temperature with respect to the renormalized elec- 
tronic level id- Both the equilibrium-state case (Vl = 0, Vr = 
0) and the steady-state case (Vl = — 1, Vk = 0) are consid- 
ered. The parameters are set as iv q — 1.5, V = 1. 



spectrum: 



A(u g ) = 



6(u> q ) 




(12) 



time (2ti/T) 

FIG. 4: (Colour online) The transient behaviours in response 
to a square-shaped pulse bias with phonon subsystem having 
a Lorentzian spectrum, (a) and (c) show Q(t) vs. t, (b) and 
(d) are Qr(t) vs. t. In (a) and (b), V = —0.5, in (c) and (d), 
Td = 3. Other parameters are the same as those in Fig. [2] 



state: 



Now the oscillations of Q(t) and Qr{t) gradually damp as 
shown in Fig. [2 Q(t) damps to a finite positive value for 
\V\ >wJor almost zero for \ V\ < uj q . Qr(t) damps along 
a slanting line for \V\ > or a horizontal line for V < 
cj°. In particular, for -u° < V < 0, Q(t) » and Q T (t) 
is negative after the system reaches a new steady state 
at the large time t [see the curves with V — —0.5 and 
V = —1.0 in Fig. [2jb)] . It means that, in the process of 
applying the sudden bias, a net energy transfer from the 
phonon subsystem to the electronic subsystem occurs. 
The phonon population n p h in the new steady state is 
less than that of the original zero-T case and it can keep 
this value for a very long time because of Q(t) w 0. To 
pay attention, the system is now in the steady state, not 
in equilibrium. So n p h being less than the original zero- 
T value does not indicate that its temperature is lower 
than 0, and this result does not breach the third law of 
thermodynamics either. 

In order to further explain and comprehend the phe- 
nomenon of heat absorption from the phonon subsystem 
at zero temperature, we have calculated the occupation 
number of zero-T phonon, n p h, and the occupation num- 
ber of electron, n e , in both equilibrium state and steady 



(n d ) = -i / G < (e)de/2TT, 



\ 2 
At: 2 



(13) 



(ei-e 2 -w 9 ) 2 uj 2 

where G<(e) iG r (e)E Q f Q f Q (e)]G a (e), 

G>(e) = -iC-(e){£ Q r Q [l _ f a ( e )]}G a (e), and 
f a (e) = l/{exp[(e - fi a )/k B T] + 1} (a = L,R). See 
Fig. [31 the occupation number of phonon behaves much 
like that of electron, which monotonically decreases 
with the elevated renormalized level of the QD, that is 
to say, n p h is approximately proportional to n c while 
changing id- At t < 0, there is no bias across the QD 
(Vl = 0, Vr = 0), so the whole system is in equilibrium. 
After applying the sudden bias, the system is driven 
out of equilibrium, until a large time, enters a new 
steady state (Vl = — 1,Vr = 0). As we can see, both 
the phonon population rtph and the electron popula- 
tion n c experience an obvious falling down from the 
equilibrium-state case to the steady-state case, it means 
that a part of zero-T phonons have indeed run away, 
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which results in the heat absorption. To be specific, for 
id = 0, the change of the average phonon number at zero 
temperature is about 0.044A^/r 2 , and the corresponding 

phonon energy is 0.066A^/r. As a matter of fact, the 
net heat absorption —Qt in real case will be a little 
less than this value, in that the electric current can also 
transfer energy to the phonon subsystem as an offset [see 
the approached level of Qt with V = — 1.0 in Fig. E^b)]. 

Finally, we deal with the square-shaped pulse bias 
Al(£), in which case Al(£) suddenly goes back to at 
t = Td- Fig. @] shows Q(t) and Qr(t) for the Lorentzian- 
spectrum phonon subsystem with different Td and V. 
While the bias going back to 0, Q(t) and Qr(t) show 
strong oscillations once again, which will gradually fade 
away. In the limit of (t —> oo), the system restores the 
equilibrium as initial, Q(t) and Qr(t) approaches zero 
and a finite positive value respectively. For \V\ > w", 
Qt(oo) is a large value that is approximately propor- 
tional to Td- But for |V| < Qt(oo) is almost zero 
regardless of Td. So if the device is operated in that range 
of | V | , there will not be much net heat generation. In par- 
ticular, as discussed above, for — cj® < V < 0, Qr{t) can 
be negative after the first pulse bias, but Qt is eventually 



> when the system recovers the equilibrium state. Ac- 
tually, because of the third law of thermodynamics, it is 
prohibited to absorb heat energy from the zero-T phonon 
subsystem under the condition of equilibrium, no matter 
how to control the bias. 

IV. CONCLUSIONS 

In summery, the transient heat generation in the QD 
driven by a step-like pulse bias is investigated. We find 
that the suddenly switched-on bias can induce a heat ab- 
sorption from the zero-temperature phonon subsystem. 
For a strong QD-lead coupling case, the current-rising 
process is much quicker than the heat generation, so the 
device can finish the electronic operation before the emer- 
gence of heat generation. In addition, for a small bias, 
the net heat generation is always very small no matter 
how to turn on and off the device. 
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